Nonlinear dispersion of stationary waves in collisionless plasmas 
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A nonlinear dispersion of a general stationary wave in collisionless plasma is obtained in a non- 
differential form from a single-particle oscillation-center Hamiltonian. For electrostatic oscillations 
in nonmagnetized plasma, considered as a paradigmatic example, the linear dielectric function is gen- 
eralized, and the trapped particle contribution to the wave frequency shift Aw is found analytically 
as a function of the wave amplitude a. Smooth distributions yield Acj ~ a^^^, as usual. However, 
beam-like distributions of trapped electrons result in different power laws, or even a logarithmic 
nonlinearity, which are derived as asymptotic limits of the same dispersion relation. 
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Introduction. — Nonlinear stationary waves, such as 
Bernstein-Green-Kruskal (BGK) modes, remain of con- 
tinuing interest [H including recently in connection 
with Raman backscattering Q and new methods of phase 
space manipulation However, essential properties of 
these waves are not apparent, because the waves are de- 
rived directly from the Vlasov-Maxwell system. The non- 
linear dispersion relations (NDR) are obtained then in a 
differential form 0, 0] , which is both specific to partic- 
ular settings and may be analytically intractable, thus 
obscuring the underlying physical picture. 

Here, we offer a universal non-differential NDR 
[Eq. ([2])] with a transparent physical meaning. The new 
NDR reveals that the nonlinear properties of a wave in 
collisionless plasma are entirely determined by one func- 
tion, namely, the single-particle oscillation-center (OC) 
Hamiltonian H [6]. Once H is found, one can study the 
nonlinear effects systematically and hence keep track of 
effects that are easy to miss in ad hoc calculations. Elec- 
trostatic waves in nonmagnetized plasma are considered 
as a paradigmatic example. For those, we show how vari- 
ous types of kinetic nonlinearities, previously known from 
different contexts, and also a new logarithmic nonlinear- 
ity are derived as asymptotic limits of the same dispersion 
relation [Eq. ([5])]. Besides that, the fundamental linear 
dielectric function is generalized [Eqs. (IT^ and (ITi|) ]. and 
the friction drag on trapped particles is predicted to af- 
fect the wave frequency sweeping in collisional plasmas. 

Basic equations. — To start, consider the plasma La- 
grangian Ls = £cm + Si^ii where Lgm is the electro- 
magnetic field Lagrangian, and Li are the Lagrangians 
of individual particles, also accounting for the interac- 
tion with the field. The plasma adiabatic dynamics on 
time scales large compared to the period of any oscilla- 
tions in the system is then governed by the time-averaged 
Lagrangian, Ls = (is)t Notice further that, in a 
stationary wave, particles can be described by some gen- 
eralized canonical coordinates Qi and momenta T'i, re- 
ferred to as OC variables [sj, such that Qi and Ti remain 
constant. Since the particle dynamics is trivial in these 
variables, let us exclude them as separate degrees of free- 
dom. This is done using Routh reduction [s, H^, which 



yields the Lagrangian L — 
the wave only. Further, since 



J2i ' 2i that describes 



(1) 



is ith particle OC Hamiltonian [9|, one obtains L — 
{Lcni)t — Hence, the wave Lagrangian L per 

unit spatial volume is given by £ = £oin — X]s"'s(^s)' 
where £om — {E"^ — B'^)x.tl{?>''^) (with averaging per- 
formed over both time and space), E and B are the 
electric and magnetic fields, summation is taken over dif- 
ferent species s, Us are the corresponding space- average 
densities, and {T-Ls) is the OC energy averaged over 'Pg- 
Assuming the wave spatial profile is prescribed, the 
dynamics of the wave is fully characterized by its am- 
plitude a (arbitrarily normalized) and canonical phase ^ ; 
by definition, the latter increases at some constant rate 
^ = oj, by 27r per the oscillation period T = juj. Hence, 
£ = £(a,^), where we used that £ cannot depend on ^ 
explicitly for it describes the dynamics on time scales 
t ^ T: cf. Rcf. 7]. In particular, varying L with respect 
to a at fixed uj yields 9a£ = [ij], or 



J_ d_ 

Stt da 



ix,t 



= 0. 
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Complemented by Eq. (P) for Ui [ij], Eq. ^ is the 
sought NDR, with advantages that it (i) applies to any 
stationary wave in collisionless plasma, (ii) has a non- 
differential form, (iii) is nonperturbative in the field am- 
plitude, (iv) allows understanding the wave properties by 
studying just %s, (v) is comprised of terms with trans- 
parent physical meaning. Below, examples are given that 
illustrate the power of this main result. 

First, revisit linear waves, in which case there clearly 
must exist modes of the form E, B oc e'*''''; then the com- 
monly known dispersion relation for linear waves without 
trapped particles should follow. To confirm this, sub- 
stitute "Hs in the dipole approximation 15[, namely, 

= -E* • a, • E/4, (3) 

where V^s^ is some function of V (and static fields, if 
any), $s is the ponderomotive potential, and cts is the 
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particle linear polarizability. Take E = ae, where e de- 
termines polarization; then daHs = — \ (e* • O-s ■ e)a and 
also i? = |n X e\a, where n = ck/w, and c is the speed of 
light. Hence, Eq. ^ gives (e* -e-e) — |n x ep = 0, where 
e = 1 + 47rns(Q!s) equals the linear dielectric tensor. 
Thus, our result coincides with the known dispersion re- 
lation at prescribed e, whereas the equation for e 
also can be recovered, by varying L with respect to e*. 

Electrostatic wave. — Now let us apply Eq. ^ to 
derive the dispersion of a nonlinear electrostatic wave 
in nonmagnetized plasma. Assume that ions are fixed; 
hence, only electron motion will be addressed, and the 
species index s is dropped. Also, neglect fluid nonlin- 
earities, which are of higher order in a than the kinetic 
nonlinearities discussed below. Then, treating the wave 
as monochromatic is anticipated to yield asymptotically 
precise description at small amplitudes [16|, ll7|. We 
hence introduce the wavenumber k and the phase velocity 
u = Lu/k. (Both u and electron velocities will be assumed 
nonrelativistic.) From Eq. ([T]), it is seen then that is 
conserved to a w-dependent term AH when transferring 
from the laboratory frame K to the reference frame K' 
where the wave field is static. Since ATi is independent 
of a, for the purpose of using Eq. ([2]) it only remains to 
find H in K', which is done as follows. 

First consider the electron true Hamiltonian in K', 



H{x,p) — /{2m) -f eipQ cos{kx), 



(4) 



where m and e are the particle mass and charge, p = 
m{v — u) is the corresponding momentum {v being the 
velocity in K) , and, for clarity, the amplitude of the po- 
tential energy is defined such that a = k"^ cLpf) / {raui'^) > 0. 
Governed by Eq. both passing and trapped parti- 
cles will undergo oscillations which are convenient to de- 
scribe in terms of the action J oc § pdx and the conjugate 
canonical phase 0, which will serve as "P and Q, in this 
case. Specifically, choose the coefficient in the expression 
for the passing-particle action such that J = \p\/k for 
large p (we assume fc > 0), and, for trapped particles, 
such that J is continuous across the separatrix. Then, 
J = Ja^/^j{r), where J — moj/k'^, and [Fig. [T](a)] 




E(r)-F(r-l)K(r), r < 1, 
i/2E(r-i), r>l 



(5) 



is a continuous function of the normalized energy r = 
{H + e(po) / i^eipo) , such that j — for a particle resting 
at the bottom of the potential trough (r = 0), with the 
corresponding value at the separatrix (r = 1) being j„ = 
A/tt. (Here K and E are the complete elliptic integrals of 
the first and second kind, respectively; cf., e.g., Ref. [ll]-) 
Since the generating function of the canonical trans- 
formation (x,p) {07 J) clearly does not depend on 
time explicitly, one gets 7^(J) = H{x,p), or 'H{JTa) — 
{2ar — a)Juj, where the dependence on a is parametric. 
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FIG. 1: Auxiliary dimensionless functions (solid): (a) j{z), 
(b) g{z), (c) ip{z), (d) q{z). The vertical dashed lines show 
where the functions are nonanalytic. Also shown are asymp- 
totes and asymptotic approximations (dashed) flowing from 
Eq. (O, except in (a), where the approximations used are 
j{z ^1)^ z + z'^/8 and j{z > 1) « 2z^^'^ - (2z^/^)"\ 



"{j) is determined by Eq. ([5]), and j 



-1/2 

defined above. Then, Eq. ^ can be rewritten as 
2 



J/J, as 
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uj =LJp■ 



g{J)F{J)dJ. 



(6) 



Here F{J) is the action distribution, and g{j) 
[dan{J,a)]/{mu^), i.e., g{j{r)) = 2r - 1 - j{r)/f {r), 



giJ) 




J < 1, 
J > 1- 



(7) 



[Notice that g{j) is continuous at the separatrix, with 
dii*) — yst with a discontinuous infinite derivative; 
Fig. IH^b).] In particular, when plasma is cold and a ~> 0, 
then all particles are passing and J = \p\/k, so 



u ± kJ/m, 



(8) 



where the sign index denotes sgn {v — u). Yet are not 
present then, and v~ <^ u, in which case Eq. ([8]) gives 
J « J [in other words, one may assume F{J) « d{J—J)]. 
Since g{j » 1) « l/(2j^), one thereby obtains g = a/2, 



meaning that Eq. ([6]) predicts uj 



as expected. 



Equation ^ describes all kinetic corrections (to the 
extent that the monochromatic-wave approximation ap- 
plies), and it readily shows how particles with given j 
affect the wave frequency. In particular, it shows that 
deeply trapped particles {j < 0.96) reduce w^, for the 
corresponding g is negative; yet those near the separatrix 
and untra ppe d ones have positive g and thus increase 
(cf. Ref. [19|). Below, we explicitly calculate for a 
number of representative cases, by formally considering 
various asymptotic expansions of the integral in Eq. ([6]) . 

Smooth distribution F{J). — First, let us assume that 
the distribution function F{J) remains finite at small J 
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or, at least, diverges less rapidly than J ^. Then, one 
can take the integral in Eq. ([6]) by parts and obtain 



2uj^ r°° 

1 1 / *(J,a)F'(J)dJ = 0, 



(9) 



where we introduced = — g{j) dJ = Ja^^'^ip{j) and 

^(j) = -J^g0)dj [Fig.lHc)], SO 



^{J, a) 




., J»Jai/2. 



(10) 



At a <C 1, 5* changes rapidly with J compared to F{J), 
if the distribution is smooth, e.g., thermal. Then, with- 
out using the explicit form of 5'( J, a) but rather drawing 
on the leading terms in Eq. (ITU|) . one can put Eq. ^ in 
the following asymptotic form: 

2 2 

e(oj,fc) + ^Cilna+-p^>*C2ai/' = 0. (11) 



Here we introduced 
e = 1 



fc6 







F'iJ)~F'iO)q 



dJ_ 

T' 



(12) 



q{z) = 1 - 2zi>{z) [Fig.lHd)], Ci = {m/kfF'{Q), C2 = 
{m/kfF"{Q), and >f = q{z) dz w 0.544. [Notice that, 
albeit determined by essentially nonlinear dynamics in 
the narrow vicinity of the resonance, q{z) nevertheless 
can affect the integrand on the thermal scale.] In partic- 
ular, when F'iQ) = 0, the nonlinear part of Eq. ([Tl]) is 
small, yielding that the nonlinear frequency shift Aw is 
also small; hence, 



(13) 



Yet, at nonzero F'{Q), the nonlinear part of Eq. pT|) 
diverges logarithmically at small a; i.e., wave interaction 
with resonant particles has a strong effect on uj. 

Equations (fTT|l - (ll3p generalize the existing NDR for 
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eigenwaves in plasmas with smooth distributions [2 
namely, as follows. First of all, notice that e, serving as a 
generalized linear dielectric function here, is a functional 
of the action distribution. Unlike the commonly used 
distribution of "unperturbed" velocities fo{v), which de- 
pends on how the wave was excited [2^, F{J) is defined 
unambiguously; thus, the above equations hold for any 
excitation scenario [while finding F{J) itself is kept as a 
separate problem]. Second, even if put in terms of fa{v), 
Eqs. pT |) -(fT3 l) cover a wider class of particle distributions. 
The latter is seen as follows. 

For example, consider a wave developed slowly from 
zero amplitude, so each J is conserved, even through 



trapping and untrapping [l^ [22, llj]- Then F{J) = 
Fo{J), index henceforth denoting the initial state. 
Yet, since there was no wave in that state, Eq. ([8]) ap- 
plies, so each ^th derivative of Fq{J) reads as Fff\j) — 
{k/mY[f^'^\v+) + {-lYf^'^\v-)]. Let us use bars to de- 
note Hmits fQ \v u±), so that fo^\v) is defined as a 
piecewise-constant function equal to the left and right 
limits for v < u and v > u correspondingly. Then, 
Ci = f^iu+) - f^iu^), C2 = f{;iu+) + J'^{u-), and 



p 



dv, 



(14) 



where q = q{\v/u — 1|). (Remarkably, contributing to e 
are both passing and trapped particles.) 

First, compare Eq. (ITil) with the usual e^ — 1 — 
i^l/k"^) - u)'^fo{v)dv, P denoting the princi- 

pal value [14| . For smooth fo{v), our e can be put in the 



same form as el, because P 



(V — uj 



'qfi{v)dv = Q. 



However, Eq. ([T4|) is valid also when /q(i') is discontin- 
uous across the resonance, a case in which el is unde- 
fined. This is because the integrand in Eq. (|14p is finite 
(piecewise-continuous), so the integral converges abso- 
lutely rather than existing only as a principal value (like 
does). Second, for smooth fo{v), when Ci = and 
C2 = 2/^'(u), Eq. (dSl) for Aw matches that in Ref. 
including the coefficient. Yet, unlike the existing theory, 
our Eqs. (fTT|) - ((T4l) apply just as well for arbitrary Ci and 
C2, in which case /o(w) may not be smooth while F{J) is. 

Beam nonlinearities. — Suppose now that, in addi- 
tion to a smooth distribution J^{J), near the resonance 
there is a phase-space clump or a hole, further termed 
uniformly as a beam with Ff){J) ^ and some average 
spatial density ni, ^ 0; namely, F{J) = J-{J) + Fh{J). 
For example, take Fi,[J) = {nh/no)S{J), where uq is the 
bulk density that enters here due to normalization. Since 
g{0) = — 1, Eq. ^ yields then, with ~ innhe^/m: 



2ujl/a, 



(15) 



[here nonlinearities due to F are neglected, and — 
uj\(ijj,k) corresponds to the linear equation], or, more 
specifically, e(cj, k) + 2uj1/{auj'^) — 0. These equa- 
tions agree with the known NDR for modes with deeply 
trapped particles [13, [HI and leads to Acj — 0(a~^), 
such that Aw < for a clump and Aw > for a hole. 

Now let nb itself depend on the wave amplitude. For 
example, a Van Kampen mode would have 2w^/(aw£) = 
rj of order one [i^; in this case, by adjusting a, any w 
can be produced for a given k, in agreement with the 
linear theory [3, 27 1. Also, consider the case when 
Ffc is constant across the trapping width: Ff,(J) = 
Ffc9(J)6( J* — J), with 9 being the Heaviside step func- 
tion. Then is proportional to the separatrix action J.^, 
[i.e., 7ib — 0(a^/^)], and one gets 



w^ = wi - [8/(37r)]a" 



^l/2„,2 



uj;jFb, 



(16) 



4 



since V'O*) = 4/(37r). Equation (IT6l) also matches the 
result found previously, e.g., in Ref. [28| . 

Finally, consider dissipation-driven effects in collisional 
plasmas. Since lo changes rapidly with small a in the 
presence of a phase-space clump or a hole, slow decay of 
a will cause frequency downshifting or upshifting, corre- 
spondingly. Yet, since the power index a in the scaling 
Aw cx a~'^ depends on how localized F},{J) is, another 
effect is anticipated, namely, as follows. Notice that a 
friction drag (say, proportional to the particle velocity) 
can cause condensation of the trapped distribution near 
the bottom of the wave potential trough [2§|. Hence, 
peaking of F{J) can occur, and a can increase gradu- 
ally up to unity. This represents a frequency sweeping 
mechanism additional to those considered in Refs. [SOj. 

Conclusions. — In summary, we show here that know- 
ing the appropriate single-particle OC Hamiltonian is suf- 
ficient to derive the fully nonlinear dispersion of a station- 
ary wave in collisionless plasma without solving Vlasov 
or Maxwell's equations. We illustrate how our theory 
reduces to results previously known from separate con- 
texts, recovering them within a single NDR. In partic- 
ular, for longitudinal electron oscillations in nonmagne- 
tized plasma, various types of kinetic nonlinearities are 
derived, including a new logarithmic nonlinearity, sim- 
ply by substituting appropriate distributions F{J) into 
Eq. Also, the linear dielectric function is generalized, 
and the friction drag on trapped particles is predicted to 
affect the wave frequency sweeping in collisional plasmas. 
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